Critical Scaling Properties at the Superfluid Transition of 4 He in Aerogel 
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We study the superfluid transition of 4 He in aerogel by Monte Carlo simulations and finite size 
scaling analysis. Aerogel is a highly porous silica glass, which we model by a diffusion limited cluster 
aggregation model. The superfluid is modeled by a three dimensional XY model, with excluded 
bonds to sites on the aerogel cluster. We obtain the correlation length exponent v = 0.73 ± 0.02, 
in reasonable agreement with experiments and with previous simulations. For the heat capacity 
exponent a, both experiments and previous simulations suggest deviations from the Josephson 
hyperscaling relation a = 2 — dv. In contrast, our Monte Carlo results support hyperscaling with 
a = — 0.2 ± 0.05. We suggest a reinterpretation of the experiments, which avoids scaling violations 
and is consistent with our simulation results. 
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Experiments on 4 He immersed in ultralight aerogel, 
forming highly porous fractal structures, have demon- 
strated new universality classes of the superfluid phase 
transition [J, Striking deviations of the critical ex- 
ponents from the bulk lambda transition appear to arise 
due to fractal correlations in the aerogel, but the details 
are not well understood. Furthermore, the new critical 
exponents tend to violate common scaling laws, such as 
the Josephson hyperscaling relation. A theoretical study 
of the unusual scaling properties of the superfluid tran- 
sition in fractal random media is motivated in order to 
better understand these surprising results. 

Experiments on helium in aerogel use ultralight aero- 
gel samples with a highly porous structure that consists 
of silica strands of about 5 nm in diameter and a dis- 
tance between neighboring strands of up to 200 nm 
The porosity can be varied in the fabrication process. 
Small-angle x-ray scattering indicates a fractal structure 
extending from a couple of nanometers up to more than 
100 nm 3]. The superfluid density of the 4 He has been 
measured by a torsional oscillator technique, and the heat 
capacity by an ac heating method Q. The superfluid 
density was found to vanish at the transition as a power 
law p s = po |*| C , where t = 1 - T/T c and C ~ 0.79 for 
m = 5% volume fraction of aerogel, C ~ 0.76 for m = 2%, 
and £ w 0.72 for m — 0.5%. These values are significantly 
greater than ( ss 0.67 found in bulk 4 He. The heat ca- 
pacity has a peak at T = T c that was fitted to the form 
C = ^4±|i|~ Q± + B, giving a + ^ a_. Alternatively, as- 
suming that a + = a- = a gives a w —0.57 for m = 5% 
aerogel, a « —0.39 for m = 2%, and a « —0.13 for 
rn = 0.5% @. Surprisingly, these exponents depend on 
the aerogel porosity, and deviate considerably from the 
value obtained from the Josephson hyperscaling relation, 
a = 2 — dv, with d = 3 and v = Q. 

The He transition in aerogel has been studied by com- 
puter simulation in Refs. [3, 1^3] , and we follow many as- 
pects of these works. Moon and Girvin Q studied a 



three dimensional (3D) XY model with aerogel modeled 
as a connected percolation fractal cluster with fractal di- 
mension Df ps 2.5. They obtained the expected scal- 
ing behavior of the superfluid density, i.e., p s ~ £ 2_d 
with d = 3, and a correlation length that diverges as 
£ ~ \T — T c \~ v with exponent v w 0.722, in agreement 
with experiments 0]. However, they did not study the 
unusual heat capacity exponents. Vasquez et al. [5J] also 
studied a 3D XY model, but with aerogel modeled with a 
diffusion limited cluster aggregation (DLCA) algorithm, 
generating a fractal cluster with dimension Df w 1.8 
[6j. They find porosity dependent exponents similar to 
those from the experiments. Specifically, the heat capac- 
ity exponent displays a similar violation of hyperscaling 
as obtained from experiments. 

In this Letter we study the superfluid transition in 
aerogel by Monte Carlo (MC) simulations and finite size 
scaling analysis, using a more extensive disorder averag- 
ing and larger system sizes than in previous simulations. 
In contrast to results from experiments and earlier simu- 
lations, we obtain critical exponents that are both poros- 
ity independent and obey usual scaling laws. In particu- 
lar, our MC data is consistent with the Josephson hyper- 
scaling relation. We reanalyze the experimental data for 
the heat capacity from Ref. in light of our new scaling 
results, and argue that the formula C = A±\t\~ a± + B 
used in the previous data analysis may not directly ap- 
ply for this problem. We present an alternative analysis 
that leads to consistency between experiments and our 
simulation results. 

DLCA model. - Following Ref. [H, we model the aero- 
gel using a DLCA algorithm 0, [i[ . We consider a simple 
cubic lattice with L x L x L sites and periodic boundary 
conditions in all directions. As a first step mL 3 "par- 
ticles" are distributed on random lattice sites, where m 
is the volume fraction of the aerogel. Particles on near- 
est neighbor lattice sites form clusters, whose structure 
is kept fixed. Clusters are allowed to diffuse by rigid 
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translations, using an iterative DLCA algorithm @. If 
two clusters meet they connect irreversibly and form a 
larger cluster. The iteration continues until a single clus- 
ter remains, which is used to model the aerogel, and plays 
the role of correlated disorder for the helium. The aero- 
gel volume fractions considered here are m — 5% and 

771 = 10%. 

3D XY model. - The universal properties of the super- 
fluid transition of 4 He are captured by the 3D XY model 
given by 



H 



(1) 



where (i, j) denotes nearest neighbor sites. The coupling 
constant Jy is set to J = 1 between pairs of 4 He sites, 
and J = if one or both lattice sites belong to the DLCA 
cluster. The excluded bonds act like quenched disorder 
for the 3D XY model, and 9i is the phase of the superfluid 
order parameter at site i occupied by helium. 

Monte Carlo simulations. - Our MC simulations use 
the collective Wolff update method lid], combined with 
the temperature exchange method [11| . which reduces 
critical slowing down of the simulation. An exchange up- 
date attempt was performed after every ten Wolff cluster 
updates. For each random bond realization we discard 
between 2 x 10 4 and 2 x 10 5 Wolff clusters to approach 
equilibrium, depending on system size, followed by the 
same number of updates collecting data. In the tempera- 
ture exchange method we simulate a set of different tem- 
peratures in parallel, and allow MC moves exchanging 
the temperatures between the different configurations. 
The results were averaged over 5 x 10 3 DLCA clusters 
for L < 50, and for 2.5 x 10 3 clusters for L > 60. We 
carefully checked for convergence of the simulation by 
increasing the number of initial discarded update steps 
until stable results were obtained. 

Calculated quantities and finite size scaling relations. 
- The superfluid density p s is proportional to the helicity 
modulus Y, which gives the increase in free energy den- 
sity as A/ = \YA 2 in the presence of a uniform phase 
twist A imposed across the system, here taken in the x 
direction [lij ■ The helicity modulus is then given by [l3| 
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where JJ = Jj+e x ,y Here (• • •) denotes thermal average, 
and [• • • ] denotes average taken over different realizations 
of DLCA clusters. 

In order to analyze MC data for finite systems of 
size L 3 , we start from the finite size scaling relation 
for the singular part of the free energy density, f s /T — 
L~ d f±(L/£), where d — 3, f± are scaling functions, and 



£ ~ |T - T c \- U is the correlation length [14J. The cor- 
responding scaling result for a phase gradient is A = 
L- 1 A±(L/£). The helicity modulus becomes Y ~ £ 2 ~ d , 
or, alternatively, Y ~ |T - T C |C [3 with ( = (d - 2)v. 
For d = 3 the finite size scaling relation becomes Q 



LY/T = (1 + aL- u ')Y{L 1/,/ t), 



(3) 



where Y is a scaling function to be determined below 
from MC data, and t = T - T c . The term aL~ u is an 
amplitude correction to scaling that will be discussed be- 
low [H. 

The Josephson hyperscaling relation for the heat ca- 
pacity exponent, a = 2 — dv 1 follows from c ~ 
d 2 f s /dT 2 ~ d 2 t dv /dt 2 ~ t dv - 2 = t~ a 14] . The heat 
capacity is given by c = [(H 2 ) - {H) 2 ] / (L 3 T 2 ), and the 
finite size scaling relation is [lq 



c = L a ' v c{L L l v t)+b, 



(4) 



where b represents the nonsingular, analytic contribution 
to the heat capacity, which depends on temperature but 
not system size. Since we focus on a narrow temperature 
interval around the transition, we take b to be constant. 
Alternatively, the scaling function can be determined di- 
rectly by eliminating b: 



(x) 



c(x,Li) - c(x,L 2 ) 



(5) 



where x is the scaling variable L x l v t and Li,L 2 denote 
two different system sizes. These two formulas gave the 
same results for the critical exponents. 

Results. - Figure [T] shows MC data for the heat capac- 
ity c (main figure) and the dimensionless helicity modulus 
LY/T (inset) vs temperature T for a set of different sys- 
tem sizes L for an aerogel volume fraction of m — 5%. 
According to the scaling relation Eq. ([3|), the dominant 
L dependence for the helicity modulus at T = T c is 
Y ~ L , showing that the superfluid transition temper- 
ature T c is where all data curves for LY/T for different 
system sizes L intersect. However, a close examination 
of the data curves for Y in Fig. [1] reveals a small sys- 
tematic drift in the intersection temperatures between 
successive pairs of system sizes, which complicates an ac- 
curate determination of T c . The drift can be included 
in the scaling analysis by assuming an amplitude correc- 
tion to scaling of the form included above in Eq. ([3|) . We 
found that varying to in the interval from 0.6 to 1 leads 
to practically indistinguishable fits to the MC data, with 
nearly the same values for the other exponents. The pre- 
cision of our MC data is not enough for a more accurate 
estimate of u>, and the results shown below are for the 
choice to = 1. For the heat capacity no correction to 
scaling was found to be necessary. 

The following method produces a finite size data col- 
lapse of MC data for the helicity modulus by fitting 
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FIG. 1: (color online) MC data for the heat capacity c (main 
figure) and the dimensionless helicity modulus LY/T (inset) 
as a function of temperature T for different system sizes L, 
with m — 5% aerogel volume fraction. 



to Eq. ([3]). The x 2 value of the fit is obtained as 
X 2 = (l/^Eoj/ dx[fi(x) - Yj{x)} 2 , where ij refer 
to n pairs of system sizes Li, Lj. The integrals are evalu- 
ated using cubic spline interpolation at 20 evenly spaced 
points within the x interval in which the pair of functions 
in the integrand overlap. The fit parameters are deter- 
mined by minimizing x 2 over a fine multidimensional grid 
of values, whose resolution is finer than the error esti- 
mates on the parameters. Data points at small values of 
Y with poor statistics are excluded from the fits. A sim- 
ilar procedure is used to fit Eq. ((4]) to the heat capacity 
data. We tried both independent scaling fits for Y and 
c, and joint fits to scale both quantities simultaneously, 
which we found to be more stable. 

To determine the critical exponents v and a we at- 
tempt some different approaches. We first assume that 
the hyperscaling relation a = 2 — dv is valid, which 
leaves the fit parameters T c , v, a, b. The result of joint 
fits of Eqs. © and flU) to MC data is shown in Fig. [2] 
for to = 5%. As seen in the figure the fit produces good 
data collapses both for Y and for c. The fit parameters 
for to = 5% are: T c = 2.1385±0.0005, v = 0.73±0.02, a = 
— 1.9,6 = 2.85, and for m = 10% (data not shown): 
T c = 2.0450 ± 0.0005, v = 0.73 ± 0.02, a = -2.0, b = 2.6. 
Error bars are estimated from the variation in the pa- 
rameter values upon varying the range of system sizes 
included in the fits. We observe a slight tendency for the 
exponents to drift towards the 3D XY results upon in- 
creasing the system sizes included in the scaling fits, but 
larger system sizes would be required to study this effect 
further. The values of T c are close to the value T c « 2.203 
for the pure 3D XY model, which is expected since the 
DLCA cluster fills only a tiny fraction of the entire system 
volume. The values of v for 5% and 10% aerogel closely 
agree with each other, and we get the corresponding heat 
capacity exponent a = 2 — iv = —0.20 ± 0.05, indepen- 
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FIG. 2: (color online) Finite size scaling data collapses of MC 
data showing the scaling functions for the helicity modulus 
Y and the heat capacity c for m = 5%. The data collapses 
giving the scaling functions are obtained from a joint fit to 
Eqs. (J3J) and Q for system sizes L > 20, assuming that a is 
given by the hyperscaling relation a — 2 — dp. 



dent of the aerogel porosity. We also tried joint scaling 
fits without assuming the hyperscaling relation for a, but 
instead treating v and a as independent fit parameters, 
with similar results. On the other hand, using the pa- 
rameter values suggested by the experiments |2j gives a 
good fit to the helicity modulus, but a very poor fit to 
the heat capacity data. Thus, from fits to our MC data 
we conclude that the exponents are porosity independent 
and obey hyperscaling, in contrast to suggestions in Rcfs. 

as. 

Comparison with experiments. - Here we compare our 
critical exponents to the experimental results in Ref . 0] • 
Our scaling of the helicity modulus gives the correlation 
length exponent v w 0.73, which compares quite well 
with the experimental values v = 0.72 — 0.79 from tor- 
sional oscillator measurements of the superfluid density. 
However our heat capacity exponent a « —0.20 is far 
from the experimental values a = —0.39 for to = 2% 
aerogel, and a = —0.57 for m — 5%. These results are 
clearly incompatible. 

We now argue that our scaling results can actually de- 
scribe the universal scaling function probed by the ex- 
periment. However, we also suggest that the experimen- 
tal data contains nonuniversal features that must be in- 
cluded in the data analyses as well. We first assume 
the presence of a finite cutoff length scale in the experi- 
mental system. Figure 3 of Ref. [2] shows that the data 
for the 4 He transition without any aerogel produces a 
sharp lambda curve. However, the presence of the aero- 
gel seems to slightly round off the maximum of the heat 
capacity curves, compared to the sharp peak predicted by 
the formula c ~ \T — T c \~ a . Rounding of the experimen- 
tal curves therefore seems to be caused by the presence 
of the aerogel, and suggests a finite length scale that may 
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FIG. 3: (color online) Experimental data of the heat capacity 
for three aerogel volume fractions, 0.5%, 2%, and 5%, repro- 
duced from Ref. [2[ (curves with symbols), together with fits 
to the scaling function (thick curves) given by Eq. JB). 



be associated with the typical aerogel pore size. 

Further, we argue that the temperature dependence of 
the superfluid condensate amplitude may contribute to 
the observed asymmetry of the experimental heat capac- 
ity curves around T c . The fit presented in Ref. [2| to 
extract a from experiments assumes that only the sin- 
gular c ~ \T — T c \~ a dependence on T is observed. We 
propose that this might be too restrictive in cases when 
a is negative, since the singular temperature dependence 
does not clearly dominate the temperature variation as 
when a > 0. We include the regular temperature depen- 
dence of the amplitude by a linear approximation of the 
form A = Aq + A\T, where Aq,Ai are constants 17 1. 
This gives 



c(T) = (A + A X T) {l*I v c\I}I v {T - T B )] + b} , (6) 

where the scaling function c is determined in Fig. [2] To 
fit Eq. ((6]) to experiments we assume a = —0.20, v = 0.73 
and use L, b, Aq, A\ as fit parameters. The results of the 
fits are shown in Fig. [3J The fit parameters are: m = 
5% : L = 835,6 = 5.03, A = 13592, A 1 = -6150; m = 
2% : L = 1107, b = 2.84, A = 65914, A 1 = -30098; m = 
0.5% : L = 1900,6 = 1.98, Aq = 193062, A x = -88428. 
We thus obtain quite good agreement in a narrow interval 
around the transition between the experiment of Ref. [2] 
and a fit to the simulated scaling function, which fulfills 
hyperscaling. The 5% data curve shows a slight deviation 
between theory and experiment on the high temperature 
side of the transition in Fig. [3] This deviation occurs 
away from the transition temperature, where additional 
temperature dependencies are expected. Note that the 
cutoff length depends roughly as L ~ 1/m 1 / 3 , which is 
compatible with the interpretation as a crossover related 
to the typical pore diameter in the aerogel. 

Discussion. - From our finite size scaling analysis of 



MC data for the 3D XY model in aerogel modeled as 
DLCA clusters, we obtain a correlation length exponent 
v w 0.73. This supports the previous conclusion 0, 0| 
that the universality class of the phase transition in the 
presence of DLCA clusters deviates from the pure 3D XY 
universality class, where v ss 0.671 [lij], even though the 
possibility of a slow crossover towards the pure 3D XY 
universality class cannot entirely be ruled out. In con- 
trast to previous results in the literature 0, 0|, our re- 
sults suggest that the critical exponents are independent 
of porosity and obey the Josephson hyperscaling relation 
a = 2 — dv, within the statistical uncertainty of the MC 
data. We argue that the formula c = A±\T-T c \~ a± +B, 
used to analyze experiments in Ref. [2| is not directly ap- 
plicable in this case. We propose an alternative analysis 
procedure based on our simulation result to fit the ex- 
perimental heat capacity data at the transition. This 
fit assumes that the experiment contains signatures of 
a finite crossover length scale, presumably related to the 
typical aerogel pore size, and of a temperature dependent 
superfluid condensate amplitude. 
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